The Color Glass Condensate and Small x 
Physics: 4 Lectures. 



Larry McLerran 

Nuclear Theory Group, Brookhaven National Laboratory, Upton, NY 11793 



Summary. The Color Glass Condensate is a state of high density gluonic mat- 
ter which controls the high energy limit of hadronic matter. These lectures begin 
with a discussion of general problems of high energy strong interactions. The infi- 
nite momentum frame description of a single hadron at very small x is developed, 
and this picture is applied to the description of ultrarelativistic nuclear collisions. 
Recent developments in the renormalization group description of the Color Glass 
Condensate are reviewed. 



1 Lecture I: General Considerations 
1.1 Introduction 

QCD is the correct theory of hadronic physics. It has been tested in var- 
ious experiments. For high energy short distance phenomena, perturbative 
QCD computations successfully confront experiment. In lattice Monte-Carlo 
computations, one gets a successful semi-quantitative description of hadronic 
spectra, and perhaps in the not too distant future one will obtain precise 
quantitative agreement. 

At present, however, all analytic computations and all precise QCD tests 
are limited to a small class of problems which correspond to short distance 
physics, or to semi-quantitative comparisons with the results of lattice gauge 
theory numerical computations. For the short distance phenomena, there is 
some characteristic energy transfer scale E, and one uses asymptotic freedom, 

a s {E) -> (1) 

as E — > oo For example, in Fig. ^, two hadrons collide to make a pair of jets. If 
the transverse momenta of the jets is large, the strong coupling strength which 
controls this production is evaluated at the px of the jet. If px » Aqcd, 
then the coupling is weak and this process can be computed in perturbation 
theory. QCD has also been extensively tested in deep inelastic scattering. In 
Fig. ^, an electron exchanges a virtual photon with a hadronic target. If the 
virtual photon momentum transfer Q is large, then one can use weak coupling 
methods. 

One question which we might ask is whether there are non-perturbative 
"simple phenomena" which arise from QCD which are worthy of further 
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Fig. 1. Hadron-hadron scattering to produce a pair of jets. 



electron 




Fig. 2. Deep inelastic scattering of an electron on a hadron. 



effort. The questions I would ask before I would become interested in under- 
standing such phenomena are 

• Is the phenomenon simple in structure? 

• Is the phenomena pervasive? 

• Is it reasonably plausible that one can understand the phenomena from 
first principles, and compute how it would appear in nature? 

I will argue that gross or typical processes in QCD, which by their very 
nature are pervasive, appear to follow simple patterns. The main content of 
this first lecture is to show some of these processes, and pose some simple 
questions about their nature which we do not yet understand. 

My goal is to convince you that much of these average phenomena of 
strong interactions at extremely high energies is controlled by a new form of 
hadronic matter, a dense condensate of gluons. This is called the Color Glass 
Condensate since 

• Color: The gluons are colored. 
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• Glass: We shall see that the fields associated with the glass evolve very 
slowly relative to natural time scales, and are disordered. This is like a 
glass which is disordered and is a liquid on long time scales but seems to 
be a solid on short time scales. 

• Condensate: There is a very high density of massless gluons. These gluons 
can be packed until their phase space density is so high that interactions 
prevent more gluon occupation. This forces at increasingly high density 
the gluons to occupy higher momenta, and the coupling becomes weak. 
The density saturates at dN/d 2 prcPrT ~ l/a^ » 1, and is a conden- 
sate. 

In these lectures, I will try to explain why the above is very plausible. 

1.2 Total Cross Sections at Asymptotic Energy 

Computing total cross sections as E — ► oo is one of the great unsolved prob- 
lems of QCD. Unlike for processes which are computed in perturbation theory, 
it is not required that any energy transfer become large as the total collision 
energy E — > oo. Computing a total cross section for hadronic scattering there- 
fore appears to be intrinsically non-perturbative. In the 60's and early 70's, 
Regge theory was extensively developed in an attempt to understand the to- 
tal cross section. The results of this analysis were to my mind inconclusive, 
and certainly can not be claimed to be a first principles understanding from 
QCD. 

The total cross section for pp and pp collisions is shown in Fig. [| 
Typically, it is assumed that the total cross section grows as ln 2 E as E — > oo. 
This is the so called Froisart bound which corresponds to the maximal growth 
allowed by unitarity of the S matrix. Is this correct? Is the coefficient of ln 2 E 
universal for all hadronic precesses? Why is the unitarity limit saturated? 
Can we understand the total cross section from first principles in QCD? 
Is it understandable in weakly coupled QCD, or is it an intrinsically non- 
perturbative phenomenon? 

1.3 How Are Particles Produced in High Energy Collisions? 

In Fig. ^ , I plot the multiplicity of produced particles in pp and in pp col- 
lisions. The last six points correspond to the pp collisions. The three upper 
points are the multiplicity in pp collisions, and the bottom three have the 
mutliplicity at zero energy subtracted. The remaining points correspond to 
pp. Notice that the pp points and those for pp with zero energy multiplicity 
subtracted fall on the same curve. The implication is that whatever is causing 
the increase in multiplicity in these collisions may be from the same mecha- 
nism. Can we compute N(E), the total multiplicity of produced particles as 
a function of energy? 
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Fig. 4. The total multiplicity in pp and pp collisions. 
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1.4 Some Useful Variables 

At this point it is useful to develop some mathematical tools. I will introduce 
kinematic variables: light cone coordinates. Let the light cone longitudinal 
momenta be 

P ± = ^(E± Pz ) (2) 

Note that the invariant dot product 

P- q=Pf qt~P + q~ -p~q + (3) 
and that 

P + p- = \{E 2 -pl) = \{p 2 T + m*)= l -ml (4) 

This equation defines the transverse mass %. (Please note that my metric 
is the negative of that conventionally used in particle physics.) 

Consider a collision in the center of mass frame as shown in Fig. [5]. In 
this figure, we have assumed that the colliding particles are large compared 
to the size of the produced particles. This is true for nuclei, or if the typical 
transverse momenta of the produced particles is large compared to Aqcd, 
since the corresponding size will be much smaller than a Fermi. We have also 
assumed that the colliding particles have an energy which is large enough so 
that they pass through one another and produce mesons in their wake. This 
is known to happen experimentally: the particles which carry the quantum 
numbers of the colliding particles typically lose only some finite fraction of 
their momenta in the collision. 




large p small p large p 

Fig. 5. A hadron-hadron collision. The produced particles are shown as red circles. 
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The right moving particle which initiates the collision shown in Fig. |^ 
has p\ <~ V2 | p z | and p± ~ iv 7 ! 771 ? 1 / \ P* \- For the colliding particles 
rriT = m pro j ec tu e , that is because the transverse momentum is zero, the 
transverse mass equals the particle mass. For particle 2, we have p\ = p\ 
andp^ =p\. 

If we define the Feynman x of a produced pion as 

x =vXlv\ (5) 

then < x < 1. (This definition agrees with Feynman's original one if the 
energy of a particle in the center of mass frame is large and the momentum is 
positive. We will use this definition as a generalization of the original one of 
Feynman since it is invariant under longitudinal Lorentz boosts.) The rapidity 
of a pion is defined to be 

V = \ln( P +/p-) = ^ln(2p+ 2 /m 2 T ) (6) 

For pions, the transverse mass includes the transverse momentum of the pion. 

The pion rapidity is always in the range —ycM < V < Ucm where 
ycM = ln(p + /mprojectue) All the pions are produced in a distribution of 
rapidities within this range. 

A distribution of produced particles in a hadronic collision is shown in 
Fig. ^. The leading particles are shown in blue and are clustered around the 
projectile and target rapidities. For example, in a heavy ion collision, this is 
where the nucleons would be. In red, the distribution of produced mesons is 
shown. 




Fig. 6. The rapidity distribution of particles produced in a hadronic collision. 



These definitions are useful, among other reasons, because of their simple 
properties under longitudinal Lorentz boosts: p^ — » K ±1 p ± where n is a 
constant. Under boosts, the rapidity just changes by a constant. 

The distribution of mesons, largely pions, shown in Fig. 4. are conve- 
niently thought about in the center of mass frame. Here we imagine the 
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positive rapidity mesons as somehow related to the right moving particle and 
the negative rapidity particles as related to the left moving particles. We de- 
fine x = P + /Pp rojecUle and x' = P~/p~ rojectae and use x for positive rapidity 
pions and x' for negative rapidity pions. 

Several theoretical issues arise in multiparticle production. Can we com- 
pute dN/dyl or even dN/dy at y = 0? How does the average transverse 
momentum of produced particles < pr > behave with energy? What is the 
ratio of produced strange/nonstrange mesons, and corresponding ratios of 
charm, top, bottom etc at y — as the center of mass energy approaches 
infinity? 

Does multiparticle production as E — > oo at y — become simple, 
understandable and computable? 

There is a remarkable feature of rapidity distributions of produced 
hadrons, which we shall refer to as Feynman scaling. If we plot rapidity 
distributions of produced hadrons at different energies, then as function of 
the distance from the fragmentation region, the rapidity distributions are to 
a good approximation independent of energy. This is illustrated in Fig. ^. 
This means that as we go to higher and higher energies, the new physics is 
associated with the additional degrees of freedom at small rapidities in the 
center of mass frame (small-x degrees of freedom). The large x degrees of 
freedom do not change much. This suggests that there may be some sort of 
renormalization group description in rapidity where the degrees of freedom 
at larger x are held fixed as we go to smaller values of x. We shall see that in 
fact these large x degrees of freedom act as sources for the small x degrees of 
freedom, and the renormalization group is generated by integrating out low 
x degrees of freedom to generate these sources. 



UN 

dy 




Fig. 7. Feynman scaling of rapidity distributions. The two different colors corre- 
spond to rapidity distributions at different energies. 
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1.5 Deep Inelastic Scattering 

In Fig. |^, deep inelastic scattering is shown. Here an electron emits a virtual 
photon which scatters from a quark in a hadron. The momentum and energy 
transfer of the electron is measured, and the results of the break up are not. 
In these lectures, I do not have sufficient time to develop the theory of deep 
inelastic scattering. Suffice it to say, that this measurement is sufficient at 
large momenta transfer Q 2 to measure the distributions of quarks in a hadron. 

To describe the quark distributions, it is convenient to work in a reference 
frame where the hadron has a large longitudinal momentum p\ adron - The 
corresponding light cone momentum of the constituent is Pjonstituenf We 
define x = ptonstituent/P hadron- This x variable is equal to the Bjorken x 
variable, which can be defined in a frame independent way. In this frame 
independent definition, x = Q 2 /2p ■ Q where p is the momentum of the 
hadronic target and Q is the momentum of the virtual photon. The cross 
section which one extracts in deep inelastic scattering can be related to the 
distributions of quarks inside a hadron, dN/dx. 

It is useful to think about the distributions as a function of rapidity. We 
define this for deep inelastic scattering as 

V = Vhadron ~ ln(l / x) (7) 

and the invariant rapidity distribution as 



dN/dy = xdN/dx 



(8) 



dN 
dy 




y 



Fig. 8. The rapidity distribution of gluons inside of a hadron. 



In Fig. |^, a typical dN/ dy distribution for constituent gluons of a hadron 
is shown. This plot is similar to the rapidity distribution of produced particles 
in deep inelastic scattering. The main difference is that we have only half of 



The Color Glass Condensate and Small x Physics: 4 Lectures. 



9 



the plot, corresponding to the left moving hadron in a collision in the center 
of mass frame. 

We shall later argue that there is in fact a relationship between the 
structure functions as measured in deep inelastic scattering and the rapidity 
distributions for particle production. We will argue that the gluon distribu- 
tion function is in fact proportional to the pion rapidity distribution. 

The small x problem is that in experiments at Hera, the rapidity dis- 
tribution function for quarks grows as the rapidity difference between the 
quark and the hadron grows. This growth appears to be more rapid than 
simply | yproj — V | or (y pr oj — y) 2 , and various theoretical models based on 
the original considerations of Lipatov and colleagues suggest it may grow as 
an exponential in | y pro j — y |-fjj (Consistency of the BFKL approach with 
the more established DGLAP evolution equations remains an outstanding 
theoretical problem.^) If the rapidity distribution grew at most as y 2 , then 
there would be no small x problem. We shall try to explain the reasons for 
this later in this lecture. 




-4 -3 -2 
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Fig. 9. The Zeus data for the gluon structure functions. 



In Fig. |9|, the Zeus data for the gluon structure function is shown. Q I 
have plotted the structure function for Q 2 = 5 GeV 2 , 20 GeV 2 and 200 GeV 2 . 
The structure function depends upon the resolution of the probe, that is Q 2 . 
Note the rise of xg(x) at small x, this is the small x problem. If one had plotted 
the total multiplicity of produced particles in pp and pp collisions on the 
same plot, one would have found rough agreement in the shape of the curves. 
Here I would use y — log(E cm /l GeV) for the pion production data. This 
is approximately the maximal value of rapidity difference between centrally 
produced pions and the projectile rapidity. The total multiplicity would be 
rescaled so that at small x, it matches the gluon structure functions. This 
demonstrates the qualitative similarity between the gluon structure function 
and the total multiplicity. 
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Fig. 10. Saturation of gluons in a hadron. A view of a hadron head on as x decreases. 



Why is the small x rise in the gluon distribution a problem? Consider 
Fig. [lOl where we view hadron head on.|||-||^] The constituents are the va- 
lence quarks, gluons and sea quarks shown as colored circles. As we add more 
and more constituents, the hadron becomes more and more crowded. If we 
were to try to measure these constituents with say an elementary photon 
probe, as we do in deep inelastic scattering, we might expect that the hadron 
would become so crowded that we could not ignore the shadowing effects of 
constituents as we make the measurement. (Shadowing means that some of 
the partons are obscured by virtue of having another parton in front of them. 
For hard spheres, for example, this would result in a decrease of the scatter- 
ing cross section relative to what is expected from incoherent independent 
scattering.) 

In fact, in deep inelastic scattering, we are measuring the cross section for 
a virtual photon 7* and a hadron, <t 7 » hadron- Making x smaller corresponds to 
increasing the energy of the interaction (at fixed Q 2 ). An exponential growth 
in the rapidity corresponds to power law growth in 1/x, which in turn implies 
power law growth with energy. This growth, if it continues forever, violates 
unitarity. The Froissart bound will allow at most ln 2 (l/x). (The Froissart 
bound is a limit on how rapidly a total cross section can rise. It follows from 
the unitarity of the scattering matrix.) 

We shall later argue that in fact the distribution functions at fixed Q 2 
do in fact saturate and cease growing so rapidly at high energy. The total 
number of gluons however demands a resolution scale, and we will see that 
the natural intrinsic scale is growing at smaller values of x, so that effectively, 
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the total number of gluons within this intrinsic scale is always increasing. The 
quantity 

,2 1 dN 

defines this intrinsic scale. Here irR 2 is the cross section for hadronic scatter- 
ing from the hadron. For a nucleus, this is well defined. For a hadron, this is 
less certain, but certainly if the wavelengths of probes are small compared to 
R, this should be well defined. If 

A 2 » A% CD (10) 

as the Hera data suggests, then we are dealing with weakly coupled QCD 
since as (A) « 1. 

Even though QCD may be weakly coupled at small x, that does not mean 
the physics is perturbative. There are many examples of nonperturbative 
physics at weak coupling. An example is instantons in electroweak theory, 
which lead to the violation of baryon number. Another example is the atomic 
physics of highly charged nuclei. The electron propagates in the background 
of a strong nuclear Coulomb field, but on the other hand, the theory is weakly 
coupled and there is a systematic weak coupling expansion which allows for 
computation of the properties of high Z (Z is the charge of the nucleus) atoms. 

We call this assortment of gluons a Color Glass Condensate. The name 
follows from the fact that the gluons are colored, and we have seen that 
they are very dense. For massless particles we expect that the high density 
limit will be a Bose condensate. The phase space density will be limited by 
repulsive gluon interactions, and be of order l/a s >> 1. The glass nature 
follows because these fields are produced by partons at higher rapidity, and 
in the center of mass frame, they are Lorentz time dilated. Therefore the 
induced fields at smaller rapidity evolve slowly compared to natural time 
scales. These fields are also disordered. These two properties are similar to 
that of a glass which is a disrodered material which is a liquid on long time 
scales and a solid on short ones. 

If the theory is local in rapidity, then the only parameter which can de- 
termine the physics at that rapidity is A 2 . Locality in rapidity means that 
there are not long range correlations in the hadronic wavefunction as a func- 
tion of rapidity. In pion production, it is known that except for overall global 
conserved quantities such as energy and total charge, such correlations are of 
short range. Note that if only A 2 determines the physics, then in an approx- 
imately scale invariant theory such as QCD, a typical transverse momentum 
of a constituent will also be of order A 2 .If A 2 >> 1/R 2 , where R is the radius 
of the hadron, then the finite size of the hadron becomes irrelevant. Therefore 
at small enough x, all hadrons become the same. The physics should only be 
controlled by A 2 . 

There should therefore be some equivalence between nuclei and say pro- 
tons. When their A 2 values are the same, their physics should be the same. 
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We can take an empirical parameterization of the gluon structure functions 
as 

1 dN A 1 '* 

where S ~ .2 — .3. This suggests that there should be the following correspon- 
dences: 

• RHIC with nuclei ~ Hera with protons 

• LHC with nuclei ~ Hera with nuclei 

Estimates of the parameter A for nuclei at RHIC energies give ~ 1 — 
2 Gev 7 and at LHC 2-3 Gev. 

Since the physics of high gluon density is weak coupling we have the 
hope that we might be able to do a first principle calculation of 

• the gluon distribution function 

• the quark and heavy quark distribution functions 

• the intrinsic pr distributions quarks and gluons 

We can also suggest a simple escape from unitarity arguments which 
suggest that the gluon distribution function must not grow at arbitrarily 
small x. The point is that at smaller x, we have larger A and correspondingly 
larger px- A typical parton added to the hadron has a size of order 1/pr- 
Therefore although we are increasing the number of gluons, we do it by 
adding in more gluons of smaller and smaller size. A probe of size resolution 
Ax > 1/pt at fixed Q will not see partons smaller than this resolution size. 
They therefore do not contribute to the fixed Q 2 cross section, and there is 
no contradiction with unitarity. 



1.6 Heavy Ion Collisions 

In Fig. O, the standard lightcone cartoon of heavy ion collisions is shown. 
To understand the figure, imagine we have two Lorentz contracted nuclei 
approaching one another at the speed of light. Since they are well localized, 
they can be thought of as sitting at x = 0, that is along the light cone, for 
t < 0. At x = 0, the nuclei collide. To analyze this problem for t > 0, it 
is convenient to introduce a time variable which is Lorentz covariant under 
longitudinal boosts 

r=y / t 2 - z 2 (12) 

and a space-time rapidity variable 
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Fig. 11. A space-time figure for ultrarelativistic heavy ion collisions. 



For free streaming particles 



, Vz , 
z — vt = —t 
E 



(14) 



we see that the space-time rapidity equals the momentum space rapidity 



If we have distributions of particles which are slowly varying in rapidity, 
it should be a good approximation to take the distributions to be rapidity 
invariant. This should be valid at very high energies in the central region. 
By the correspondence between space-time and momentum space rapidity, 
it is plausible therefore to assume that distributions are independent of r\. 
Therefore distributions are the same on lines of constant r, which is as shown 
in Fig. [n]. At z = 0, t = t, so that r is a longitudinally Lorentz invariant 
time variable. 

We expect that at very late times, we have a free streaming gas of 
hadrons. These are the hadrons which eventually arrive at our detector. At 
some earlier time, these particles decouple from a dense gas of strongly in- 
teracting hadrons. As we proceed earlier in time, at some time there is a 
transition between a gas of hadrons and a plasma of quarks and gluons. This 
may be through a first order phase transition where the system might exist 
in a mixed phase for some length of time, or perhaps there is a continuous 
change in the properties of the system 

At some earlier time, the quarks and gluons of the quark-gluon plasma 
are formed. This is at RHIC energies, a time of the order of a Fermi, perhaps 
as small as .1 Fermi. As they form, the particles scatter from one another, 
and this can be described using the methods of transport theory. At some 
later time they have thermalized, and the system can be approximately de- 
scribed using the methods of perfect fluid hydrodynamics. 



V = V 



(15) 
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In the time between that for which the quarks and gluons have been 
formed and r = 0, the particles are being formed. This is where the initial 
conditions for a hydrodynamic description are made. 

In various levels of sophistication, one can compute the properties of 
matter made in heavy ion collisions at times later than the formation time. 
The problems are understood in principle for r > Tf orma ti on if perhaps not 
in fact. Very little is known about the initial conditions. 

In principal, understanding the initial conditions should be the simplest 
part of the problem. At the initial time, the degrees of freedom are most 
energetic and therefore one has the best chance to understand them using 
weak coupling methods in QCD. 

There are two separate classes of problems one has to understand for 
the initial conditions. First the two nuclei which are colliding are in single 
quantum mechanical states. Therefore for some early time, the degrees of 
freedom must be quantum mechanical. This means that 

AzAp z > 1 (16) 

Therefore classical transport theory cannot describe the particle down to 
r = since classical transport theory assumes we know a distribution function 
/(p, x, t), which is a simultaneous function of momenta and coordinates. This 
can also be understood as a consequence of entropy. An initial quantum state 
has zero entropy. Once one describes things by classical distribution functions, 
entropy has been produced. Where did it come from? 

Another problem which must be understood is classical charge coher- 
ence. At very early time, we have a tremendously large number of particles 
packed into a longitudinal size scale of less than a fermi. This is due to the 
Lorentz contraction of the nuclei. We know that the particles cannot interact 
incoherently. For example, if we measure the field due to two opposite charge 
at a distance scale r large compared to their separation, we know the field 
falls as 1/r 2 , not 1/r. On the other hand, in cascade theory, interactions are 
taken into account by cross sections which involve matrix elements squared. 
There is no room for classical charge coherence. 

There are a whole variety of problems one can address in heavy ion 
collisions such 

• What is the equation of state of strongly interacting matter? 

• Is there a first order QCD phase transition? 

These issues and others would take us beyond the scope of these lectures. 
The issues which I would like to address are related to the determination 
of the initial conditions, a problem which can hopefully be addressed using 
weak coupling methods in QCD. 

1.7 Universality 

There are two separate formulations of universality which are important in 
understanding small x physics. 
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The first is a weak universality. This is the statement that physics should 
only depend upon the variable |7| 

As discussed above, this universality has immediate experimental 
consequences which can be directly tested. 

The second is a strong universality which is meant in a statistical me- 
chanical sense. At first sight it appears to be a formal idea with little relation 
to experiment. If it is however true, its consequences are very powerful and 
far reaching. What we shall mean by strong universality is that the effective 
action which describes small x distribution function is critical and at a fixed 
point of some renormalization group. This means that the behavior of corre- 
lation functions is given by universal critical exponents, and these universal 
critical exponents depend only on general properties of the theory such as 
the symmetries and dimensionality. 

Since the correlation functions determine the physics, this statement says 
that the physics is not determined by the details of the interactions, only by 
very general properties of the underlying theory! 



2 Lecture II: A Very High Energy Nucleus 

In this lecture, I will consider the properties of a single nucleus. I will 

develop the theory of the small x part of the nucleus, the components most 
relevant in the high energy limit. I will begin with some general consider- 
ations. This will largely be done to develop approximations which will be 
useful later, and leads directly to the Color Glass description. I then present 
a brief review of light cone quantization. Finally, I turn to a computation of 
the color fields which describe the nuclear wavefunction at small x. I show 
that in a simple approximation for the Color Glass, one recovers saturation, 
and that the phase space density of the fields is of order 1 / a s which is typical 
of a Bosc condensate. 



2.1 Approximations and the Color Glass 

In the previous lecture, I argued that when we go to small x 

A * = ^R^» A %cn (18) 

that the theory is weakly coupled a s (A) << 1. The typical transverse mo- 
mentum scale of constituents of this low x part of the hadron wavefunction 
is 

pl~ A 2 » 1/R 2 had (19) 
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This equation means that the scale of transverse variation of the hadron 
is over much larger sizes than the transverse De Broglie wavelength. I can 
therefore treat the hadron as having a well defined size and collisions will 
have well defined impact parameter. 

For our purposes, it is sufficient to treat the hadron as a thin sheet of 
infinite transverse extent. The transverse variation in radius can be reinserted 
in an almost trivial generalization of these considerations. The thinness of the 
sheet follows because I shall assume that the sources for the fields at small 
x come from partons at much larger x which are Lorentz contracted to size 
scales much smaller than can be resolved. In Fig. [l|, a nucleus in the infinite 









Fig. 12. A single nucleus in the infinite momentum frame as seen by a small x 
probe. 



momentum frame is shown, within the approximations described above. 

Recall from the first lecture, we introduced rapidities associated with 
produced particles in hadron-hadron collisions, 

y = Un{p + /p-) = ln(V2p + /M T ) 
= ln(V2p+ ad /M T ) + Hp + /p+ ad ) ~ y had - ln(l/x) (20) 

This expression shows that the rapidity of produced hadrons can be written 
in the form used to describe the rapidity of constituents of a hadron. If we 
were to think of both the constituents and produced partons as pions, they 
would be the same, or alternatively if we think of both the produced and 
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constituent partons as gluons. We can convert to spacetime rapidity using 
the uncertainty relation p x^ ~ 1. 

y ~ ^ln(x'/x + ) ~ y had - ln{x~pl ad ) (21) 

We have assumed in deriving this relationship that the typical values of the 
proper time r = V 'x+x~ are not large compared to natural scales such as 
a transverse mass. These relations argue that all rapidities, up to shifts of 
order one, are the same. We can identify all momentum-space and space- 
time rapidities! This has the profound consequence that at high energies 
momentum space and space-time are intrinsically correlated, and particles 
which arise from a localized region of momentum space rapidity also arise 
from a localized region of space-time rapidity. 

Now we illustrate a high energy hadron in terms of space-time rapidity. 
This has the effect of spreading out the thin sheet shown in Fig. |l2], as 
shown in Fig. [l^. Note that the partons which are shown have an ordering 
in momentum space rapidity which corresponds to their coordinate space 
rapidity. Fast partons are to the left. In this Figure, I have drawn a tube of 




^max y 

Fig. 13. A single nucleus shown in terms of the space-time rapidity. The red circles 
indicate partons. 



transverse extent dx which goes through the nucleus. I take dx << 1 Fm so 
that one is resolving the constituents of ordinary hadrons. Notice that when 
dx — > 0, the longitudinal separation between hadrons which intersect the 
tube becomes large. If I also require that the energy is high enough so that 
there is always a large number of partons which intersect the box (which are 
longitudinally well separated) , then one can think of a source associated with 
the charge inside the box, and this charge has a random distribution over the 
transverse area of the box. (At what scale the source becomes random is not 
entirely clear from this discussion. This will be resolved more carefully later.) 

In the limit where l/yl<<da;<<l Fm, there are many charges inside 
the box of dimension dx. The charge should go over to a classical charge on 
this resolution scale because we can ignore commutators of charges 

I [Q a ,Q b ] |=| if abc Q c \« Q 2 (22) 
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We can define a current associated with this charge which is localized in the 
sheet as 

JZ = 5» + 5(x-)q a (x T ) (23) 

The + component of the current is the only important one because the sheet 
is traveling near the speed of light. The source q a (xr) is a c-number color 
charge density which is a random variable on the sheet. It is only defined on 
scales 1/A « dx « 1 Fm. The delta function of x~ expresses the fact that 
the source is on a thin sheet. In fact, for many applications, we will have to 
relax the delta function assumption, and work with a charge density which 
includes the effect of distribution in x~ as 

q a (x T ) = J dx- p a (x-,x T ) (24) 

and where for many purposes 

p a (x-,x T )~6(x-)q a (x T ) (25) 

We now know how to write down a theory. It is a theory where one 
computes the classical gluon field in stationary phase approximation and 
then integrates over a random source function. Its measure is 

Z = J[dA][dp]exp^iS[A] +iJ+A- J dx~ d 2 x T p2{ ^^^ (26) 

In this theory, we have assumed that the sources are randomly distributed as 
a Gaussian. This turns out to be an approximation valid in a particular range 
of resolution dx, and can be fixed up for a wider range. This will be discussed 
when we do the renormalization group. The sources and fields are coupled 
together in the standard J • A form. This results in the problem that the 
extended current conservation law D^J^ = makes J not an independent 
function. This problem can be avoided by introducing a generalization of 
the J • A coupling. This generalization turns out to be important for the 
renormalization group analysis of this theory, but is not important when we 
compute the classical field associated with these sources. 

The above theory implicitly has cutoffs in it. We have discussed the 
range of dx for which this effective theory is valid. Implicit in the analysis is 
that the fields we are computing have p + values much less than those of the 
sources. This implies there is an upper p + cutoff in the fields A considered. If 
we were to Fourier analyze the sources p, they would have their support for 
| p + | which is greater than that of the cutoff. This cutoff is of course entirely 
arbitrary, and the lack of dependence of physical quantities upon this cutoff 
forms the basis of the renormalization group. 

Notice that this theory, in spite of having a gauge dependent source, is 
gauge invariant on account of the integration over all sources. This computa- 
tion of classical fields associated with sources and then averaging over sources 
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is similar to the mathematics of glasses. The physical origin of this similarity 
is the Lorentz time dilation of the source for the fields and the disorder of 
the gluon field. The Lorentz time dilation is of course an approximation, and 
if one were to observe these classical fields over long enough time scales they 
would evolve, as do the atoms in a glass. 
Notice that 

< p a {x) P \y) >= S ab S^(x - y)p 2 (x-) (27) 

so that ji 2 is the charge squared per unit transverse area per unit x~ scaled 
by l/(iV c 2 -l). 



2.2 Light Cone Quantization 

Before discussing the properties of classical fields associated with these 
sources, it is useful to review some properties of light cone quantization. [fiof 
This will allow us to pick out physical observables, such as the gluon density, 
from expectation values of gluon field operators. 
Light cone coordinates are 

x ± = -L(x°±x 3 ) (28) 



and momenta 



P 



± 



1 



(P°±P 3 ) (29) 



V2 V 

The invariant dot product is 

p ■ x — p t ■ x t — p + x~ — p~x + (30) 

where p t and xt are transverse coordinates. This implies that in this basis 
the metric is <? H = g h = —1, g l i = where i,j refer to transverse 
coordinates. All other elements of the metric vanish. 

An advantage of light cone coordinates is that if we do a Lorentz boost 
along the longitudinal direction with Lorentz gamma factor 7 = cosh(y) then 
p± _> e ±y p ± 

If we let x + be a time variable, we see that the variable p~ is to be inter- 
preted as an energy. Therefore when we have a field theory, the component 
of the momentum operator P~~ will be interpreted as the Hamiltonian. The 
remaining variables are to be thought of as momenta and spatial coordinates. 
In Fig. |TJ, there is a plot of the z,t plane. The line x + = provides a sur- 
face where initial data might be specified. Time evolution is in the direction 
normal to this surface. 

We see that an elementary wave equation 

(p 2 + M 2 )(f> = (31) 
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is particularly simple in light cone gauge. Since p 2 — p 2 — 2p + p this equation 
is of the form 

« = ^ <W 

is first order in time. In light cone coordinates, the dynamics looks similar to 
that of the Schrodinger equation. The initial data to be specified is only the 
value of the field on the initial surface. 

In the conventional treatment of the Klein-Gordon field, one must specify 
the field and its first derivative (the momentum) on the initial surface. In light 
cone coordinates, the field is sufficient and the field momentum is redundant. 
This means that the field momentum will not commute with the field on the 
initial time surface! 

Lets us work all this out with the example of the Klein Gordon field. 
The action for this theory is 

1 



\ (33) 



The field momentum is 

n {xt , X -) = ^=8^=^ (34) 

Note that 27 is a derivative of <p on the initial time surface. It is therefore 
not an independent variable, as would be the case in the standard canonical 
quantization of the scalar field. 

We postulate the equal time commutation relation 

[n(x t ,x-),4>{iH,y-)] = ~S^ 3 \x - y) (35) 
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(The factor of two in the above expression is subtle and comes from a careful 
reduction of constrained Dirac bracket quantization for the classical theory to 
quantum field theory. It can be checked by verifying that we get the correct 
result for the Hamiltonian.) Here the time is x + = y + = in both the the 
field and field momentum. We see therefore that 

d-[<j>{x)^{y)} = - l -5VXx-y) (36) 

or 

[<Kx),4>(y)] = - l -e(x~ y-)5W{x y) (37) 

Here e(v) is 1/2 for v > and -1/2 for v < 0. 

These commutation relations may be realized by the field 

d 3 p 



f 



-{e lpx a{p)+e- lpx a)(p)} (38) 



/ p +>o (2^) 3 2vV 
Using 

[a(p),aHq)}=2p+(2n) :i S^(p- q ) (39) 

one can verify that the equal time commutation relations for the field are 
satisfied. 

The quantity l/p + in the expression for the field in terms of creation 
and annihilation operators is singular when p + = 0. When we use a princi- 
ple value prescription, we reproduce the form of the commutation relations 
postulated above with the factor of e(x~ — y~). Different prescriptions corre- 
spond to different choices for the inversion of -J=r . One possible prescription is 
the Leibbrandt-Mandclstam prescription l/p + = p~ /(p + p~ + ie). This pre- 
scription has some advantages relative to the principle value prescription in 
that it maintains causality at intermediate stages of computations and the 
principle value prescription does not. In the end, for physical quantities, the 
choice of prescription cannot result in different results. Of course, in some 
schemes the computations may become prohibitively difficult. 

The light cone Hamiltonian is 

f d3 P Pt+M 2 u . , . fAn . 

with obvious physical interpretation. 

In a general interacting theory, the Hamiltonian will of course be more 
complicated. The representation for the fields in terms of creation and anni- 
hilation operators will be the same as above. Note that all particles created 
by a creation operator have positive P + . Therefore, since the vacuum has 
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P + = 0, there can be no particle content to the vacuum. It is a trivial state. 
Of course this must be wrong since the physical vacuum must contain con- 
densates such as the one responsible for chiral symmetry restoration. It can 
be shown that such non-perturbative condensates arise in the P + = modes 
of the theory. We have not been careful in treating such modes. For pertur- 
bation theory presumably to all orders, the above treatment is sufficient for 
our purposes. 



2.3 Light Cone Gauge QCD 

In QCD we have a vector field A£. This can be decomposed into longitudinal 
and transverse parts as 

At = ±={A° a ±Al) (41) 

and the transverse as lying in the two dimensional plane orthogonal to the 
beam z axis. Light cone gauge is 

A+ = (42) 

In this gauge, the equation of motion 

D^F"" = (43) 
is for the + component 

D,F t+ - D+F-+ = (44) 
which allows one to compute A~ in terms of A 1 as 

A- = ^D l d+A l (45) 

This equation says that we can express the longitudinal field entirely in terms 
of the transverse degrees of freedom which are specified by the transverse 
fields entirely and explicitly. These degrees of freedom correspond to the two 
polarization states of the gluons. 
We therefore have 

K(x) = J {2 ^ 3 P 2p+ (e ipx 4(p) + e-**a?(p)) (46) 
where 

K(P), <Ho)] = 2p + S ab ^(27rfS^(p - q) (47) 
where the commutator is at equal light cone time x + . 
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2.4 Distribution Functions 

We would like to explore some hadronic properties using light cone field 
operators. For example, suppose we have a hadron and ask what is the gluon 
content of that hadron. Then we would compute 

d ^^= <h \a)[p)a{p)\h> (48) 

If we express this in terms of the gluon field, we find 

d -^r = ^r <h 1 ^<p.* + )^°(-p.* + ) i h > ( 49 ) 

which can be related to the gluon propagator. The quark distribution for 
quarks of flavor i (for the sum of quarks and antiquarks) would be given in 
terms of creation and annihilation operators for quarks as 

^=<h\ {b\(p)bi(p) + dj(p)d t (p)} | h > (50) 

where b corresponds to quarks and d to antiquarks. The creation and annihi- 
lation operators for quarks and gluons can be related to the quark coordinate 
space field operators by techniques similar to those above. pT| 



2.5 The Classical Gluon Field 

To compute the gluon distribution function, we need the expectation value 
of the gluon field. To lowest order in weak coupling, this is given by com- 
puting the classical gluon field and then averaging over sources. The classical 
equation of motion is 

= 5 v +p{x-,x T ) (51) 

To solve this equation we shall work in the gauge A~ = 0, and then 
gauge rotate the solution back to lightcone gauge A + = 0. The solution in 
A~ = gauge is 

A 1 = 

-V|^ + = p (52) 

Here ~p = W(x)pU(x) is the source which has been gauge rotated to this new 
gauge. Since the measure for integration over sources is gauge invariant, we 
do not have to distinguish between these sources since we can rotate one into 
the other. 

To rotate back to lightcone gauge we use 



9 



(53) 
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so that the gauge rotation matrix U is 

d+U = -igUA + (54) 
where 

A + =a = -^Vp (55) 
The solution is||-@ 

= Pexplig I dz~ a(z~ , xt) { (56) 



There is a choice of boundary condition here associated with Xq . The am- 
biguity with this choice is associated with a residual gauge freedom. We 
shall resolve this by choosing retarded boundary conditions, x^ — > — oo. This 
boundary condition lets us construct the solution for U at some knowing 
only information about a for x~ < Xi . 

The solution in light cone gauge is therefore 

A+ = A- = 

A 1 = -C/V'C/t (57) 
9 

If x~ is outside of the range of support of the source p, this can be written 
as 

A^Wi-i-yv'yt (58) 

9 

where 

/POO \ 

dz~ a(z~ , xt) ) (59) 



V'(xt) = Pexp ( ig I 
\ J — t 



We now have an explicit expression for the gluon field in terms of the 
sources .[p|,[^2[ For our Gaussian weight function, we can now compute the ex- 
pectation value of the gluon fields which gives the gluon distribution function. 



The details of such a computation are given in Ref. 12 . It is a straightforward 
computation to perform: One can expand out the exponentials and compute 
term by term in the expansion. The series exponentiates. One subtlety occurs 
due to logarithmic infrared infinity which is regulated on a scale of order of 
a Fermi where transverse charge correlations go to zero since all hadrons are 
color singlet . [ji~3| . The result is 

N 2 — '\ 1 

< A> a (x,x + )A* a (0 7 x + ) > = ^r-j x 



Tra s N c xip 



[1 - exp{x 2 T Q 2 Mx 2 T A 2 QCD )/±\) (60) 
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In this equation, the saturation momentum is defined as 

Ql = 2jrN c a 2 s J dx-^ 2 (61) 

and is of the order a 2 times the charge squared per unit area. 

This expression can be Fourier transformed to produce the gluon dis- 
tribution function, with result as shown in Fig. EM We can understand this 




plot from the properties of the coordinate space distribution function. We 
notice that the dominant scale factor in the problem is Q s , so to a first ap- 
proximation everything scales in terms of this quantity. Large pr corresponds 
to small xt, and the coordinate space distribution behaves as ln(x T ) which 
corresponds to l/p T . This is typical of a bremstrahlung spectrum. At larger 
xt, distribution is of order l/x T , which Fourier transforms into ln{p T ) at 
small pt- The softer xt dependence at large xt can be traced to a dipole 
cancelation of the fields. The monopole charge field, seen at short distances 
is ln(x T ) and the dipole cancelation should set in at large distances when 
one cannot resolve individual charges, and reduce this by two powers of Xt- 

The overall scale of the curve is l/a s . The quantity we are plotting is 
in fact the phase space density of gluons. At small a s , this density becomes 
large, and the Color Glass becomes a condensate. Hence the name, Color 
Glass Condensate. 

This form of the gluon distribution function illustrates how the problems 
with unitarity can be solved. Let us assume that the saturation momentum 
Q s is rapidly increasing as x — > 0. If we start with an x so that Q » Q s , 
then as x decreases, the number of gluons which can be seen in scattering 
rises like Q 2 . (See Fig. |l^.) Eventually, Q s becomes larger than Q, in which 
case, the number of gluons rise slowly, like ln(Q s ). At this point the cross 
section saturates since the number of gluons which can be resolved stops 
growing, and we are consistent with unitarity constraints. 
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The gluon distribution function is defined to be 



xG(x,Q 2 ) = fPTjr^T 
J d 2 p T dy 



(62) 



This behaves in the saturation region as irR 2 Q 2 , and in the large Q region as 
irR 2 Q 2 . We expect that Q 2 S ~ charge 2 / area and due to the random nature of 
the way charges add, Q 2 ~ R. Therefore in the saturation region, the gluon 
distribution function is proportional to the surface area of the hadron, that 
is the gluons can only be seen which are on the surface of the hadron. In the 
large Q region, one sees gluons from the entire hadron, that is, the hadron 
has become transparent. 

2.6 The Structure of the Gluon Field 

The gluon field arises from a charge density which is essentially delta func- 
tion in x~ . In order to solve the equations of motion, the field must have a 
discontinuity at x~ = 0. This can be achieved with a field which is a two 
dimensional gauge transform of zero field strength on one side of the sheet 
and a different gauge transform of zero on the other side. The field strength 
is therefore zero if jU and v are both in the two dimensional transverse 
space. If cither index is — , it also vanishes since there is no change in the x + 
direction. The only non- vanishing component is therefore F l+ , and this is a 
delta function in the x~~ direction. Since F 4± = E % ± e^fP, we see that 



The fields are therefore transversely polarized to the direction of motion 
and live in the two dimensional sheet where the charges sit. These are the 
non-abelian generalizations of the Lienard-Wiechert potentials of electrody- 
namics. The density of these fields is of order l/a s . A picture of the Color 
Glass Condensate is shown in Fig. 

3 Lecture III: Hadron-Hadron Collisions and the 
Initial Conditions for Heavy Ion Collisions 

3.1 Phenomenology of Mini-Jets 

In the last lecture, we argued that at small x, the typical gluon constituent 
of a hadron acquires a transverse momentum of order Q s , and that this can 
grow as x — > 0. This leads us to hope that in hadron-hadron collisions, this 
will be the typical momentum scale of particle production. If true, then the 
processes are weakly coupled and computable using weak coupling methods. 



(63) 
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This is reminiscent of past attempt to compute particle production by 



mini-jets 1 14 On dimensional grounds, the cross section for jet produc- 
tion da l 'dycPpr ~ a l/Pr- If we attempt to compute the total cross section 
for jet production, 



da 2 f d Pt 

dy ~ as J A » QOD Pt 



the result is infrared sensitive, and presumably would be cutoff at Aqcd- In 
early computations, one introduced an ad hoc cutoff which was fixed, and 
hopefully large enough so that one could compute the minijet component. 
This of course left unanswered many questions about the origin of this cutoff, 
and the effects of particles produced below the cutoff scale. 

In this lecture, we will argue that the Color Glass Condensate cuts off 
the integral at a scale of order Q s , the saturation momentum. At large pt, 
dimensional arguments tell us that the density of produced particles has the 
form 



1 dN 1 



K - 



irR 2 d 2 pTdy a s p\ 



(65) 



The factor of l/p T comes about because the high pr tail is controlled by per- 
turbation theory. The 1 / a s arises because of the large density of gluons in the 
condensate. In fact, if we can successfully formulate the particle production 
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problem classically, we expect that in general 



1 



dN 



ttR 2 cPpxdy 



-f{QI/p 2 t) 



(66) 



At large px » Q s , F ~ Q a s /Pt an< ^ f° r smau Pt « Qs, F should be slowly 
varying (logarithmic) or a constant. A plot is shown in Fig. [l7| 



7tR 2 dy 




A QCD Q sat 



Fig. 17. The pr distribution for mini-jets produced by a Color Glass Condensate. 



A word of caution should be injected about the interpretation of mini- 
jet production. Typically it is assumed that there is a simple relationship 
between the multiplicity of produced gluon jets and the multiplicity of pions. 
Usually, Npi on is taken to be some constant of order one times N g i uon . In 
our considerations, we can only talk about the gluon mini-jet production, 
and it is beyond the scope of these lectures to relate this to the final state 
multiplicity. Suffice it to say that the situation is controversial, particularly 
in heavy ion collisions where there can be much final state interaction. jl^] . 

Recall that in heavy ion collisions, we expect that Q 2 ~ A 1 / 3 . At large 
pr, Eqn. pa predicts that 



dN . . ^Qi . . Ai/3 

d 2 prdy p% p? 



vR 2 ^ „ " (67 ) 



This result is consistent with hard incoherent scattering. At small pt, 

"' A ttR 2 ~ A 2 / 3 (68) 



d 2 p T dy 



which is consistent with much shadowing, and the gluons are produced from 
the surface of the nuclei. 

The total multiplicity per unit rapidity 



dN o2 f ,2 Q 



^4 



If / d' PT ^ ~ R'Qi ~ A (69) 
dy J Q 2 p\ 
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is proportional to A, just as in color string models! This is because for the 
Color Glass Condensate, the cutoff in transverse momentum depends on A. 
(If one was careful with the factors of a s in the above equation, one would 
predict mild logarithmic modifications of the linear dependence on A.) In 
addition to the A dependence, there is also a correlation between the energy 
dependence of the gluon distribution function at saturation and the multi- 
plicity of minijets since 

ttR 2 Q 2 s = [ dx'G{x',Q s ) (70) 



a relationship which follows from the last lecture. 

We can be a little more careful with the numerical factors which deter- 
mine the saturation momentum. (Tt]] Using the results of last lecture, 

2 2irN c a 2 s 2 

Qs = ^2(^2 _ color ['*■) 

Here Q 2 co i or is the color charge squared of all quarks and gluons at larger x 
values than that of interest. For a quark, 

Qluark = Nc{N 2_ 1) tV T a = ^ 

and for a gluon 

Q gluon — (*^y2 _ ^) 

(73) 

Wc find that 

✓~)2 ^ quark . ^c-^ gluon /7/1\ 

Q color - + (7V 2 ~ \ <V 

If we plug in numbers, at RHIC energies corresponding to x ~ 1CP 2 , Q s ~ 
1 - 2 GeV 



3.2 Classical Description of Hadron Collisions 

We want to describe the collision of two ultra-relativistic hadrons. A collision 
is shown in Fig. [l8|. The hadrons have been Lorentz contracted to thin sheets 
and a Color Glass Condensate sits in the planes of both sheets. 

Before the collision the non-zero fields are for right moving nucleus, 
F' l+ ~ S(x~) and for the left moving nucleus F l ~ ~ 5(x + ). Before the nuclei 
pass through one another, nothing happens and the fields in each sheet are 
static. When they pass through one another, the sum of these two fields is 
not a solution of the equations of motion, unlike the case in electrodynamics, 
and this induces a time evolution of the fields. |l8[ 

One can understand this from the vector potentials. In Fig. O a space- 
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Fig. 18. A collision of two ultra-relativistic hadrons. 




Fig. 19. A space-time diagram for the vector potentials in hadron-hadron scatter- 



time diagram is shown for the scattering. In the backward light cone, Region I, 
the field is a pure two dimensional gauge transform of zero field. In crossing 
into Regions II and III, the fields must have a discontinuity to match the 
charge on the surfaces of the lightconc. This requires the vector potential 
to be different gauge transforms of zero field strength, G2 and G3 in these 
regions. Now in going to Region IV, one could solve either for the sources on 
the left edge of the forward light cone with a gauge transform of zero or the 
right edge of the forward light cone with a different gauge transform of zero. 
One cannot satisfy the equations of motion for the fields in the presence of 
the sources on both edges of the light cones with the same gauge transform of 
zero field strength. One must produce a field in the forward light cone which 
is not a gauge, and therefore matter is produced. 

The situation in QCD is completely different than that in electrodynam- 
ics. In electrodynamics, one must produce pairs of charged particles to make 
matter in the forward light cone. This arises from a quantum correction to 
the equations of motion. In QCD, matter is produced classically. 

The procedure for solving this problem is now straightforward, in prin- 
ciple. One solves the classical equations of motion in the forward light cone 
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with boundary conditions at the edges of the forward light cone. We will 
shortly determine these boundary conditions and the form of the solution in 
the forward light cone. Then one evolves the equations of motion into the 
far future. At some time, the energy density becomes dilute, and the field 
equations should linearize in some gauge. One can then identify the quanta 
of the linearized fields in the standard way that one does classical radiation 
theory in electrodynamics. 

3.3 The Form of the Classical Field 

Before the collision, the form of the classical field can be taken as 
A+ = A- = 

A i = 9(x-)8(-x + )a\(x T ) + 6(-x-)6(x+)ai(x T ) (75) 

where the a 1 are two dimensional gauge transforms of zero field. We will 
consider the collision of identical hadrons. The solution in the forward light 
cone is therefore expected to be boost invariant. After the collision, a boost 
invariant solution is 

A + = x + a(r, Xt) 
A~ = x~/3(t, xt) 

A l = al(T,x T ) (76) 
We can choose the gauge 

x+A- + x-A+ = (77) 
so that 

a(T,x T ) = -/3(T,XT) (78) 
In the forward light cone, the equations of motion are 

-^d T r 3 d T a - [D\ [D\ a}} = (79) 

and 

-d T rd T ai - igr 2 [a, [D\a\] - [D j ,F ji ] = (80) 

T 

The boundary conditions are determined by matching the solution in 
Regions II and III to that in the forward light cone. The result is that a and 
ol\ must both be regular as r — > and 

aKO, x T ) = a\(x T ) +a l 2 (x T ) 

a(0,x T ) = ^-[a\(x T ),ai(x T )] (81) 
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The problem is now well denned, and these equations may be numerically 
solved. 

The behaviour of these solution at large r can be extracted. With V(xt) 
an element of the group, the solution is a small fluctuation field up to a 
possible large gauge transformation 

oi(t,Xt) = Ve(r, xt)V^ 

ai(T,x T ) = V(ei(T,x T )+ l -d i )Vi (82) 
The small fluctuation fields e and e % solve the equations 

-^d T T 3 d T e - V^e = (83) 

and 

-d T rd T e l - {V 2 T 5 ij - V'V J V = (84) 

T 

At large r, these linear equations can be Fourier analyzed with the result 
and 

In these equations ui ~\ k T \. 

One can compute the energy distribution associated with these fields as 

5b-w?" , " )P (87) 

and the multiplicity distribution is given by dividing this by u, that is 

^ = SFp*" ! < 88 » 

These last two formulae correspond to those of free quantum filed theory when 
we replace a(p), a*(p) by the creation and annihilation operators a(p), a'(p). 
The a(p) are the classical quantities which correspond to the quantum cre- 
ation and annihilation operators. These formulae show how to use the clas- 
sical solutions to compute distributions of produced minijets. 
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3.4 Numerical Results for Mini-Jet Production 

Krasnitz and Venugopalan have numerically solved the classical equations for 
mini-jet production. jl9] This involves finding a gauge invariant discretization 
of the classical equations of motion. One then solves the classical equations 
for a fixed p\ and p2, and extracts the produced radiation. An ensemble of 
sources are produced with the Gaussian weight of the Color Glass, which 
then produces an ensemble of radiation fields. These fields are then averaged 
to generate the mini-jet distributions. 

In Fig. the form of the numerical results for mini-jet production is 
illustrated. At large px, the results of analytic studies are reproduced which 
up to logarithms is ~ 1/Px- At px < Q s , the distribution flattens out.[i~g||- 
| jl7| To good numerical accuracy, the result in this region can be fit to a two 



dN 



\ I" (P t IQs) 




QCD 



Fig. 20. An illustration of the results generated by numerical simulation of the 
classical equations for mini-jet production. 



dimensional Bose Einstein distribution, 
1 dN 1 1 



irR 2 d 2 p T dy a s e^vrlQs - 1 



(89) 



where k is a constant of order 1. 

The result at large px can be computed analytically by expanding the 
equations in powers of the gluon field. At high px, the phase space is not 
so heavily occupied, so a field strength expansion makes sense. At small px, 
it is not at all certain that the result is in fact an exact two dimensional 
Bose- Einstein distribution . Ipjf - fHf In any case, the origins of this distribution 
have nothing at all to do with thermodynamics, and it is a useful example 
of the traps one can fall into if one assumes that exponential distribution 
corresponds to a temperature and thermalization. 



34 Larry McLerran 



3.5 pA Scattering 

An interesting example of minijet production is given by the collision of 
two hadrons of different size.||2^|-P^| We will generically refer to this as pA 
scattering, although most of our considerations could be generalized to A' A 
nuclear collisions. In Fig. ^l], the transverse momentum distribution for pA 




Fig. 21. The pr distribution for particles produced in a pA collision. 



scattering is shown. 

There are three distinct regions, which follow from the fact that there 
are two saturation scales, Qf and Q p , and Q p « Qf since (Qf) 2 ~ A 1 / 3 . 
At very large px where the fields from both nuclei are small, the distribution 
can be computed from perturbation theory, and the distribution falls as 1 /pj, 
and is proportional to (QfQ p ) 2 . This first region is for pt » Qf- An 
intermediate region where the field from the nucleus is strong but the field 
from the proton is weak and can be treated perturbatively. This intermediate 
region is for Q p « pr « Qf- There is finally the region where px « Q p s 
where both fields are strong. 

We expect that in the intermediate region, the transverse momentum 
dependence will be in between the flat behaviour at small px and the 1 jp\ 
behaviour characteristic of large pt, The naive expectation is l/pf* in the in- 
termediate region. The total multiplicity can be computed if one understands 
this intermediate region since the dominant contribution arises here. The 
strength in this intermediate region should involve the total charge squared 
from the proton, but that from the nucleus should go like p\ so that when 
combined with the \jp\, one gets a distribution proportional to \jp\. This 
softer behaviour of the distribution function follows since we are inside the 
region where we expect coherence from the field of the nucleus, and since the 
distribution should extrapolate between 1 /pj, at very large pt and a constant, 
up to logarithms, at very small pt- 

In fact, it is possible to compute the behaviour in this intermediate 
region. The equations for the classical production can be analytically solved 
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for any pt » Qg- The solution in the forward light cone are plane waves 
which are gauge transformed by the field of the large nucleus. The boundary 
conditions determine the strength of these waves. 

For the total multiplicity, in the large px region dN '/ 'dycPpr ~ A 1 / 3 . We 
expect that as we interpolate between the proton fragmentation region and 
that of the nucleus, we go between O(l) and 0(A 1 ^ 3 ) as shown in Fig. p2[ 
For pr in the intermediate region, we expect that dN/dyd 2 pT is of order 1 



Fig. 22. The distribution in rapidity for dN/dyd 2 pr in a pA collision. 



except for a small region of rapidity around the fragmentation region of the 
nucleus. The total integrated multiplicity arises from this latter region so we 
expect that dN/dy - O(l). 

3.6 Thermalization 

After the gluons are produced in hadron-hadron collisions, they may rescatter 
from one another. If one goes to very small x so that the density of gluons 
becomes very large, one expects that the gluons will eventually thermalize. 
Due to the very large typical pr, ot s << 1, arid this takes a time r ~ \/(a 2 Q s ) 
which is longer by a factor of 1/a 2 than the natural time scale. The system 
therefore becomes dilute relative to its natural scale. 

In the first diagram of Fig. there is ordinary Coulomb scattering. 
When all processes which populate and depopulate phase space are summed, 
this diagram is only naively logarithmically divergent, and is cutoff by the 
density dependent Debye mass. p g i uons « p\- In the second diagram, there 
is no such cancellation, and the diagram is of order l/(a syj /p g i uons ). At a 
time of order l/(a 2 s Q s ) for a density decreasing like 1/r as we expect for 
ultrarelativistic nuclear collisions, the diagram is enhanced by a factor of 
l/a s . This cancels the extra factor of a s coming from the diagram being 
higher order in perturbation theory. |^| . 

What appears to happen is that as the system gets more dilute, it ther- 
malizes due to multigluon production. This will modify the relationship be- 
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Fig. 23. The diagrams for gluon scattering which lead to thermalization. 



tween the number of gluons produced as mini-jets and the pion multiplicity. 
How this actually works is not yet fully understood. 

4 Lecture 4: The Renormalization Group 

The effective action for the theory we have described must be gauge invariant 
and properly describe the dynamics in the presence of external sources. For 
the theory which we have written down in past lectures with the J- A coupling 
of source to field, gauge invariance is only retained if we impose 

D„ = (90) 

This equation presents a problem in our formulation since it implies that the 
source cannot be independently specified from the field. This did not present a 
problem for the classical theory since one could find a solution which solved 
the constraint. When we compute quantum corrections and proceed to a 
renormalization group treatment, we must be more careful. 
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In a clever series of papers, p5[ it was shown that one can generalize the 
J • A coupling. This lead Jalilian-Marian, Kovner, Leonidov and Weigert to 
propose the action 

S = -\f d 4 xF* u Fr + ±J d 2 x t dx-5{x-) 

xp a (x t )trT a exp\i [ dx + T ■ A-(x)\ (91) 



In this equation, the matrix T is in the adjoint representation of the gauge 
group. This is required for reality of the action. When this action is extremized 
to get the Yang-Mills equations, one can identify the current and show that 
the current is covariantly conserved. This action is invariant under gauge 
transformations which are identified at x + = ±oo. (Even this can be corrected 
to get a fully invariant theory if one generalizes even further to complex time 
Keldysh contours. As shown in Ref. p7| , this further generalization does 
not affect the renormalization group in lowest non-trivial order.) This is a 
consequence of the gauge invariance of the measure of integration over the 
sources p. This will be taken as a boundary condition on the theory. In general 
if we had not integrated over sources, one could not define a gauge invariant 
theory with a source, as gauge rotations would change the definition of the 
source. Here because the source is integrated over in a gauge invariant way, 
the problem does not arise. 

In the most general gauge invariant theory which we can write down is 
generated from 

Z = [ [dp]e- F ^ f [dA]e iS t A >d (92) 



This is a generalization of the Gaussian ansatz described in the previous 
lecture. It allows for a slightly more complicated structure of stochastic vari- 
ation of the sources. The Gaussian ansatz can be shown to be valid when 
evaluating structure functions at large transverse momenta. 

F G aussian[p\ = ^ J dx~d 2 X t ~^T~\ ( 93 ) 

This theory is an effective theory valid only in a limited range of rapidity 
much less than the rapidity of the source. The sources for this theory sit at 
higher rapidity. This happens because as we go to lower values of rapidity, the 
fluctuations in the field are integrated out and are replaced by sources and an 
integral over fluctuations in the source. The renormalization group equations 
which we will describe are what make the theory independent of this cutoff. 
To fully determine F in the above equation demands a full solution of these 
renormalization group equations. This has yet to be done, although there are 
now approximate solutions for small and large transverse momentum of the 
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We can understand this a little better by imagining what happens when 
we compute a quantum correction to the classical theory. This quantum cor- 
rection will generate terms proportional to a s /n(yl + /p + ) where A + is the p + 
cutoff for our effective theory. Clearly these corrections are small and sensi- 
ble only if e~ 1 /° lB A + << p + << A + . If we want to generate a good effective 
theory at smaller values of p + , we need to break the theory into intervals of 
p + with each interval sufficiently small so that the quantum fluctuations are 
small and computable. The relation between one interval and the next is the 
renormalization group. 

The remarkable thing that happens when one integrates out the fluctu- 
ations interval is that only the function F which controls the source strength 
is modified! The functional form of F is modified so that this equation is of 
the form 

^Z = -H{p,5/8p)Z (94) 
dy 

where 

y = ln(A+/A+) (95) 
and Af f are the cutoffs at the initial and final values, and 

Z = e- F (96) 

This equation is of the form of the time evolution for a two spatial dimension 
quantum field theory where the coordinates are p and the momenta are d/ dp. 



4.1 How to Compute the RG Effective Hamiltonian 

In the Eqn. ^M], the renormalization group Hamiltonian H was introduced. 
I will here outline how it is computed. We first take the theory defined for 
p + < Aj . We integrate out the quantum fluctuations. In particular, the two 
point function is 

g*(x,y) =< {A\x) + 8A\x)){A\y)+8A\y)) > (97) 

In this equation, A is the classical background field and SA is the small 
fluctuation. At the momentum scales which will be of interest for p + < A~f, 
it is sufficient to consider the equal time limit of this correlation function. 
We now identify 

< SASA > = G< Sp5p > G 

= G X G (98) 

where G is the Greens function in the classical background field A. We also 
identify 

< 8 A > = G < Sp > 

= a (99) 
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We can get exactly the same result by modilying the weight function so that 
we reproduce x an d <J and move the cutoff to Aj so that there are no longer 
quantum fluctuations to integrate out. This is the origin of the form of Eqn. 

E 

Some technical comments about the computations are required. One 
must be extremely careful of gauge. The gauge prescriptions of retarded or 
advanced for l/k + singularities are used. We were not able to effectively use 
either Leibbrandt-Mandelstam or principle value prescriptions although this 
may be possible in principle. When one computes propagators in background 
fields, one gets analytic expressions in terms of line ordered phases of the 
source p. It is most convenient to compute these in SA~ — gauge and 
express things in terms of the source in A~ = gauge, and then rotate 
results back to lightcone gauge. This can be carefully done only when the 
l/k + singularity is properly regularized. 

If we change variable to space-time rapidity, we can define 

1 

a{y,x T ) = - 

and 



a(y, x T ) = —^zpiy, x T ) (100) 



V 2 T 



V\y,x T )=Pexp(ig j" dy'a{y' 7 x T ^j (101) 
After much work, one finds 

H = y J d 2 x T J la (x T )J la {x T ) (102) 
where 

J ia (x T )= f ^^4(l-1/t(^ r)% , T ))^ 6 (103) 
J tt [x — z) z ioa D (y,ZT) 

The Hamiltonian is positive definite and looks like a pure kinetic energy term 
(up to the multiplicative non-linearities) with no potential. p7| 

The renormalization group above can also be seen to be a consequence of 
equations for correlation functions of V(y, x T ).^^ In Ref. |2jJ it was shown 
that these equations for correlation functions were almost the same as those of 
the renormalization group. There was an error in this analysis associated with 
the subtleties of gauge fixing, and when repaired gives that these equations 
are precisely equivalent. ^). Meanwhile, Weigert showed that the equations 
for the correlation functions could be summarized as a Hamiltonian equation 
of the form above. |3(| which was also shown to be precisely the equations for 
the renormalization group Hamiltonian. p7[ 



4.2 Quantum Diffusion 



The Hamiltonian presented in the previous section is analogous to that with- 
out a potential. If we were to ignore the non-linearities associated with the 
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matrices V, this would be the Hamiltonian for a free theory with only mo- 
menta and no potential. 

If there was a potential in the Hamiltonian, then at large times, the 
solution of the above renormalization group equations would be trivial, 

Z ~ exp{-yE ) (104) 

where E a is the ground state energy. All expectation values would become 
rapidity independent and the solution to the small x problem would be trivial: 
x independence. 

The solution to the above equation is more complicated. One can see 
this by studying a one dimensional quantum mechanics problem: 

d v 2 

±Z = -*-Z (105) 
with solution 

Z = -^expi-x 2 /2y) (106) 

This equation describes diffusion. The width of the Gaussian in x grows with 
time. This is unlike solving 

*.Z = -(£ + V{x))z (107) 

In this latter case, the coordinate x evolves towards the minimum of V, and 
then does undergo small fluctuations around this minimum. 

We see therefore that the non-triviality of the small x problem in QCD 
arises because of the quantum diffusive nature of the renormalization group 
equations. 



4.3 Some Generic Features of the Renormalization Group 
Equation 

If we compute the correlation function of two sources using Eqn. |94|, we find 
that 

-f < p(x)p(y) >=-< p(x)p(y)H > (108) 

ay 

At large kx when the fields are linear, the gluon structure function is the 
same as the source-source correlation function up to a trivial factor of l/k T . 
(The momenta kx is conjugate to the coordinate xt — Ut-) If we ignore the 
non-linearities in H , keeping the lowest order non- vanishing terms, and if we 
integrate by parts the factors of 5/5a(y,XT), we get a closed linear equation 
for the correlation function. This is precisely the BFKL equation. 

In fact, in the region where the equations are linear, one is in the high kx 
limit, and this also reduces to the DGLAP and BFKL equations, which are 
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known to be equivalent if one computes distribution functions to leading order 
in a s ln{l/ x)ln(Q 2 ), where Q is some typical momentum for the correlation 
function, Q ~ kx- When the non-linearities are important, the non-linearities 
of this equation cannot be ignored. 

The situation is as shown in Fig. G4. In the linear region, one can choose 



Color Glass Condensate 



ln(l/x) 



ln(A QCD ) 

111(0*2) 

Fig. 24. The various regions of evolution for structure functions in the ln(l/x)- 
ln(Q 2 ) plane. 



to evolve using linear equations. In the ln(Q 2 ) direction, the equation is 
the DGLAP equation and in the ln(\/x) direction, it is the BFKL equation. 
There is a boundary region in the ln(l/x)-ln(Q ) plane. Within this boundary 
region, there is a high density of glue and the evolution becomes non-linear. 
One always collides with this region if one decreases x and holds Q 2 fixed or 
decreases Q 2 holding x fixed. 

4.4 Some Limiting Solutions of the Renormalization Group 
Equations 

In the small kx region, we expect that correlation functions such as 
< V(x)V'(y) > are very small, since we are probing the theory at distance 
scales long compared to natural correlation lengths. In this limit, one might 
be able to ignore the non-linearities in the renormalization group equations. 
Using that x l /x 2 — V'/V 2 , we have then 

/<iV J^ T) ~ JM, < , | -ij | > j^yj^, (109) 

The solution for Eqn. [m| is 

F=^- dyd 2 XTV T a{y,XT)V T a{y,x T ) (HO) 
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The small hr functional F is a pure scale-invariant Gaussian. It is universal 
and independent of initial conditions. 

In the large px region, we perform a mean held analysis. The result is 
that discussed in Lecture II. 

For details of the analysis leading to the results of this section, the in- 
terested reader is referred to 



4.5 Some Speculative Remarks 

The form of the renormalization group equation appears to be simple. It 
looks like it might even be possible to find exact solutions. In remarkable 
works, (28), |3l) Balitsky and Kovchegov have shown that the equation for the 
two line correlation function W(x)W^(y) where W is in the fundamental 
representation becomes a closed non-linear equation in the large N co i or limit. 
This means that at large N co i or one can compute this correlation function at 
arbitrarily small x including all the non-linearities associated with small x. 

Although Kovchegov's original derivation was for large nuclei, the result 
can be shown to follow directly from the renormalization group Eqn. This 
is done by taking the expectation value of < W(x)W'(y) >, using the form 
of the Hamiltonian and a factorization property of expectation values true in 
large N co i or . A derivation is presented in Ref. ^(| for the interested reader. 

This result is interesting in itself since it means that all of the saturation 
effects for i*2(x, Q 2 ) may be computed at small x. The Balitsky-Kovchegov 
equation has been solved numerically. p2| 

More important, it suggests that perhaps, at least in large N co i or , the 
full renormalization group equations may be solved for F. 



5 Concluding Comments 

In these lectures, constraints of space and time have forced me to not men- 
tion many of the exciting areas that are currently under study. One of these 
areas, is diffraction. [^3)-|3j]. One can show that the same formalism which 
gives deep inelastic scattering also gives diffraction and that there is a simple 
relation between diffractive structure functions and deep inelastic scatter- 
ing. I have also not developed a formal treatment of deep inelastic scattering 
within the Color Glass Condensate picture. [j35|-|3r|. 

The last lecture is very sketchy, and should provide an introduction to 
the literature on this problem. The derivation of the results discussed in that 
lecture are onerous, and all the details have been omitted in these lectures. 
In some sense this is good, since the most interesting part of this problem 
is to understand and solve the renormalization group equations, and at least 
this problem is clearly stated, and free from the technical details from which 
it arises.. 
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An area which should be better understood from the perspective pre- 
sented above is the nature of shadowing for nuclei at small x. This relates 
deep inelastic scattering and diffraction in a non-trivial way, and the Color 
Glass Condensate is one of the few theories available which pretends to treat 
both consistently. 

The other area where there is much potential is the production of quarks 
in hadron-hadron collisions. In particular, the charm quark may provide us a 
real clue about non trivial dynamics since its mass is very close to the scale 
of the Color Glass Condensate for large nuclei at accessible energies. 
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